Astronomy & Astrophysics manuscript no. roos'arXiv 


© ESO 2008 


February 4, 2008 





A model of accelerating dark energy in decelerating gravity 

M. Roos 

Departments of Physical Sciences and Astronomy, FIN-00014 UNIVERSITY OF HELSINKI 
e-mail: matts dot roos at Helsinki dot Fi 

Preprint online version: February 4, 2008 

ABSTRACT 

Context. The expansion of the Universe is accelerated as testified by observations of supernovae of type la at varying redshifts. 
Explanations of this acceleration are of two kinds: modifications of Einstein gravity or new forms of energy, coined dark energy. An 
example of modified gravity is the braneworld Dvali-Gabadadze-Porrati (DGP) model, an example of dark energy is Chaplygin gas. 
Both are characterized by a cross-over length scale r c which marks the transition between physics occurring on our four-dimensional 
brane, and in a five-dimensional bulk space. 

Aims. Assuming that the scales r c in the two models are the same, we study Chaplygin gas dark energy in both self-accelerating and 
self-decelerating flat DGP geometries. The self-accelerating branch does not give a viable model, it causes too much acceleration. 
Methods. We derive the Hubble function and the luminosity distance for the self-decelerating branch, and then fit a compilation of 
192 SNela magnitudes and redshifts. This determines a confidence region in the space of the three parameters of the model. 
Results. Our model with the self-decelerating branch in flat space fits the supernova data as successfully as does the ACDM model, 
and with only one additional parameter. 

Conclusions. In contrast to the ACDM model, this model needs no fine-tuning, and it can explain the coincidence problem. It is 
unique in the sense that it cannot be reduced to a cosmological constant model in any other limit of the parameter space than in the 
distant future. If later tests with other cosmological data are successful, we have here a first indication that we live in a five-dimensional 
braneworld. 
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1. Introduction 

The demonstration by SNela that the Universe is undergoing an 
accelerated expansion has stimulated a vigorous search of mod- 
els to explain this unexpected fact. Since the dynamics of the 
Universe is conventionally described by the Friedmann equa- 
tions which follow from the Einstein equation in four dimen- 
sions, all modifications ultimately affect the Einstein equation. 

The left-hand-side of the Einstein equation encodes the geo- 
metry of the Universe in the Einstein tensor G^, the right-hand- 
side encodes the energy content in the stress-energy tensor T MV . 
Thus modifications to G^y imply some alternative geometry, 
modifications in T MV involve new forms of energy densities that 
have not been observed, and which therefore are called dark en- 
ergy. 

A well-studied model of modified gravity is the Dvali- 
Gabadadze-Porrati (DGP) braneworld model (Dvali |2000l 
Deffayet 2001) in which our four-dimensional world is a FRW 
brane embedded in a five-dimensional Minkowski bulk. The 
model is characterized by a cross-over length scale r c such that 
gravity is a four-dimensional theory at scales a <K r c where mat- 
ter behaves as pressureless dust. In the self-accelerating DGP 
branch gravity "leaks out" into the bulk at scales a » r c and 
the cosmology approaches the behavior of a cosmological con- 
stant. To explain the accelerated expansion which is of recent 
date (z ss 0.5 or a « 2/3), r c must be of the order of 1. In the 
self-decelerating DGP branch gravity "leaks in" from the bulk at 
scales a » r c , counteracting the observed dark energy accelera- 
tion. 

Another well-studied model introduces into the density 
p v and pressure p v of a fluid called Chaplygin gas (Kamenshchik 



2001, Bilic 2002) following historical work in aerodynamics 
(Chaplygin I1904I I. This model is similar to the DGP model in 
the sense that it is also characterized by a cross-over length scale 
below which the gas behaves as pressureless dust, and above 
which it approaches the behavior of a cosmological constant. 
This length scale is expected to be of the same order of magni- 
tude as the r c scale in the DGP model. 

Both the self-accelerating DGP model in flat space and the 
standard Chaplygin gas model have problems fitting present su- 
pernova data, as demonstrated by Davis (120071 1. In the standard 
Chaplygin gas model the Jeans instability of perturbations be- 
haves like CDM fluctuations in the dust-dominated stage {a <K 
r c ), but disappears in the acceleration stage (a » r c ). The com- 
bined effect of suppression of perturbations and non-zero Jeans 
length leads to a strong ISW effect and thus of loss of power 
in CMB anisotropics (Amendola 2003] Bento [2003). This has 
led to generalizations to higher-dimensional braneworld models 
which appear less motivated, and which require more parame- 
ters. 

We, instead, combine the standard DGP model with the stan- 
dard Chaplygin gas model. 

This paper is organized as follows. In Section 2 we dis- 
cuss how to identify the cross-over scales in the DGP and 
Chaplygin gas models. This idea is motivated by the similar- 
ities in the asymptotic properties of the models, and was first 
presented in (Roos 2007). In Section 3 we discuss the flat-space 
self-accelerating basic DGP model with and without standard 
Chaplygin gas dark energy, and in Section 4 we turn to the self- 
decelerating DGP model combined with standard Chaplygin gas. 
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In Section 5 we summarize our results, and in Section 6 we dis- 
cuss them and conclude. 

2. Cross-over scales 

On the four-dimensional brane in the DGP model, the action of 
gravity is proportional to M 2 p( whereas in the bulk it is propor- 
tional to the corresponding quantity in 5 dimensions, Ml. The 
cross-over length scale is defined as 



3. DGP gravities with and without Chaplygin gas 

Let us now return to the Friedmann equation ([3]) and solve it for 
the expansion history H(a). Substituting Q. r from Eq. Q , p v (a) 
from Eq. d9J, and Q,„ = Q^aT 3 , it becomes 



H(a) 



Cl n + Q° cT 3 + Q A ^1 + (AD. r a 2 ) 



-3 



1/2 



■(10) 



M 2 
2M 3 



(1) 



Note that Q, v and do not evolve with a. In the limit a <K r c 
this equation reduces to two terms which evolve as gT 3 , thus 
behaving as dust with density parameter QjJ, + Q / i/(4Q ri ,) 3 . In 
the limit a » r c , Eq. ([T0| describes a model with a cosmological 



It is customary to associate a density parameter with this, 

= (2r c H r 2 - (2) 

The Friedmann equation in the DGP model may be written 
(Deffavet EOOTt 



constant Qa = — JO, rc + -y/O lv + Q A - 

At present, when a = 1 and H = Ho, we solve it for Q^, 



= 1 - 2e - Jl + (4£\.a 2 ) ~ 3 . 



(11) 



// 2 



1 



// 2 



/cp 



(3) 



In the well-studied standard self-accelerating DGP model, e 
+ 1 and D. A - 0, so that 



87rG/3, and p is the total cosmic Q° =1-2 JCl r .. 



where a = (1 + z) , , 
fluid energy density p = p m + p v . Clearly the standard FRW 
cosmology is recovered in the limit r c — > oo and p^ — > 0. In the 
following we shall only consider A: = flat geometry. The self- 
accelerating branch corresponds to e = +1, the self-decelerating 
branch to e = - 1 . 

Since ordinary matter does not interact with Chaplygin gas, 
one has separate continuity equations for the energy densities p m 
and p^ respectively. In DGP geometry the continuity equations 
have the same form as in FRW geometry (Deffayet 2001), 

p + 3H(p +p) = 0. (4) 

Pressureless dust with p — then evolves as p m (a) oc gt 3 . 

The Chaplygin gas pressure is p v = -A/p v , where A is a con- 
stant with the dimensions of energy density squared. The conti- 
nuity equation for Chaplygin gas is then 



(12) 



py + 3H p^ 

\ Pip 

which integrates to 



0, 



(a) = VA + B/c 



(5) 



(6) 



P<p 

where B is an integration constant. Thus this model has two free 
parameters. Obviously its limiting behavior is 



VB /BV'" IB\ 

p^a) oc — for a <K I — 1 , p v (a) oc -p for a » I — I .(7) 

Om first central assumption is that the DGP model cross- 
over scale r c and the Chaplygin gas cross-over scale (B/A) 1 ^ 6 
are about the same. If we choose the proportionality 

(B\ U6 r— 

(-) =r c H =2^K c , (8) 

this permits rewriting Eq. (|6]l in a dimensionless form. Making 
use of Eq. Q, p^ becomes 



<B\ 1 ^ 6 



Pip (a) = hIk- 1 D. a ^1 + (4D. r a 2 ) 



(9) 



where we have replaced the energy density *\[A by the dimen- 
sionless density parameter = Hq 2 k VA. 

The identification of the two cross-over scales evidently re- 
duces the number of free parameters in Eq. ([3]) by one: they are 
Q. rc , Q.a and Q. m — Kp m H^ 2 . 



This equation represents the condition for flatness, and corre- 
sponds to the linear relation QjJ, + Qa = 1 in the KCDM model. 
Here, however, Eq. ( 12 1 is nonlinear, causeing Qj), always to be 
smaller than in the ACDM model. This is the reason why the 
standard self-accelerating DGP model is a worse fit to SNela 
data th an the ACDM model [cf. eg. Davis (120071) . Rydbeck 
(I2007H . The failure has led to studies of various generalized 
DGP models implying higher-dimensional bulk spaces and ad- 
ditional free parameters that detract from its original simplicity 
and elegance. 

The inclusion of the Chaplygin gas term with > in 
Eq. (Hi leads a further reduction in the value of £2", and thus 
to an even worse fit to SNela data. DGP self-acceleration and 
Chaplygin gas dark energy simply yield too much acceleration, 
separately as well as in combination. In the next Section we 
therefore turn to what represents our second central assumption, 
self-decelerating DGP gravity with e - -1. Its expnasion history 
is still given by Eq. ( 10 1 and its flat-space condition by Eq. ( flT) . 
The physics then changes, as one can see best in Eq. (Ill, where 
the two last terms get opposite signs. (The case with Q. A = is 
not interesting here, because it does not lead to any acceleration.) 



4. Data and method of analysis 

The data we use to test this model are the same 192 SNela as 
in the compilation used by Davis (2007 ) which is a combination 
of the "passed" set in Table 9 of Wood-Vasey (2007) and the 
"Gold" set in Table 6 of Riess d2007b . 

We are sceptical about using CMB and BAO power spectra, 
because they have been derived in FRW geometry. SNela data 
are, however, robust for our analysis, since the distance moduli 
are derived from light curve shapes and fluxes, that do not de- 
pend on the choice of cosmological models. In one of our fits we 
nevertheless include a value for Q° as a Gaussian prior which 
we take from Table 2 of Tegmark (2006 ), who has obtained it in 
a multi-parameter fit to WMAP and SDSS LRG data. Tegmark's 
value is Q° = 0.239 + 0.018/ - 0.017, but we do not use these 
lcr errors which have been obtained by marginalizing over all 
other parameters, and which would constrain our fit too strongly. 
We take the Q. m prior to have a large error, AO" = 0.09, in or- 
der not to bias the conclusions from the SNela data set. We do 
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not marginalize, but quote full three-dimensional confidence re- 
gions: a lcr error then corresponds to a contour & x\ est + 3.54 
around the best value v? ,. 

^ best 

The Davis' compilation lists magnitudes pi, magnitude er- 
rors AyU, for SNela at redshifts Zi, i = 1, 192. We compute model 
magnitudes 



fi(zi, ci m , o rc , Q. A ) = 5 Log[d L (z ; , n m , n Ic , n A )] + 25 , 

where the luminosity distance in Mpc at redshift z, is 

dz 



d L (zi,Sl m ,n r c ,R A ) = (1 +z, 



Jo 



H(z) 



(13) 



(14) 



-1. 



We then search in the parameter space for a minimum of the 



where H(z) is given by Eq. ( |T0| > with e 

We 
^ 2 sum 

A ' 2 = g| /i, '" //(z/ ' nB "" fc ' n>t) ) 2 1 ( °- 24 - Q '" ] 2 (15) 



Occasionally we do not include the last term. 

The calculations are done with the classical CERN program 
MINUIT (James and Roos ll975t which delivers x\ s v P arameter 
errors, error contours and parameter correlations. 



5. Results 

Our first fit to determine the region in O^, £l r( . , O^-space, re- 
quired by the SNela data, keeps all the parameters free, except 
that Qj), is restricted to positive values and an upper limit is im- 
posed on Q A . We do not include the last term in Eq. ( 15 I. We 
find a solution at 



OL = 0.40 



+0.13 

at lower limit 



, Q,, = 132$*, Q A = 3.7 a X Pper " mit (16) 



with x 2 = 195.1. We have checked that the standard AC DM 
model fits the same data with the (insignificantly) higher value 
X 1 = 195.6 (as was also found by Davis 2007J- 

Since the confidence region exceeds the limits fixed for Qj)i 
and £l A , obviously y 1 is ver y insensitive to values of these pa- 
rameters near their limits. Moreover, the |(Q° , Q r .)| correlation 
coefficient is 0.97, thus one of those parameters is almost su- 
perfluous. In that sense, our model appears almost as a two- 
parameter model. One could actually fix at an arbitrary value, 
but that would be an assumption ad hoc. 

To cure these fitting problems, we proceed instead to in- 
clude the last term in Eq. ( 15 i as a weak prior, choosing AD. m 
as large as possible, in this case AQ", = 0.09 . This has the ef- 
fect of neatly putting the error contours well inside all imposed 
(and now unnecessary) limits, and reducing the correlation co- 
efficient, |(Q°,Q r< )| to 0.87. We now find as best solution the 
parameter values 



Cl° m = 0.26 ± 0.16, R n = 0.82^, Q A = 2.21+°f> 



(17) 



with x 2 = 195.5. 

In Fig.l we plot the best fit confidence region in the 
(ilJJ,, D, v )-plane, a banana-shaped closed contour . A cross in the 
Figure marks the point of best fit. 



6. Discussion and conclusions 

To learn how well our three-parameter model fits the confidence 
region determined by the data, we turn to the flat-space condi- 
tion, Eq. (fTTl). Here this condition is a surface in the fijj,, Q r ,, Q A - 
space which, if our model is successful, should cut the banana- 
shaped confidence region in Fig.l. Unfortunately the exact value 
of Q. A is not known, so we must draw Eq. (Hi for several values. 
Obviously the model is a good fit to the observational data when 
Qa is within the lcr range quoted in Eq. 



17) 



Substituting the solution from Eq. ( 17 1 into Eq. |2]), the value 
of the cross-over scale is r c = 0.55/ y/Ho- The relation Eq. ^ 
was a conjecture that could well have been different by some 
numerical factor. Then the fitted values of the free parameters 
would have changed, but a good fit could still have been ob- 
tained. Thus we do not consider that the relation Eq. |8]l is a 
fine-tuning condition. In contrast to the ACDM model, to the 
Quintessence model, and to many other models, the present 
model does not imply any fine-tuning. 

To compare models by using Information Criteria BIC 
(Schwarz [19781 or AIC (Akaike [F774b as do Davis (120071 we 
consider extremely crude. The reason is that no information on 
parameter correlations is included. If one parameter pair has a 
correlation coefficient near 0.99, it should be counted as a single 
parameter. 

It is easy to explain the coincidence problem in this model. 
It is caused merely by the ratio of the scales of the action, the 
Planck scale M P i on our brane and the bulk scale M5. These 
constants just happen to have particular time-independent values 
which determine the DGP cross-over scale r c by the definition 

Our model should still be tested against other cosmological 
data, but this has to wait until CMB and BAO power spectra have 
been derived for five-dimensional braneworld cosmology. As for 
ISW data, the problems encountered in the simplest Chaplygin 
gas model are alleviated if not eliminated by the presence of 
DGP self-deceleration. 

If this model meets all criteria, we have here a first indication 
that we live in a five-dimensional braneworld. 

Note added in print. Recently cosmologies embedding the 
generalized Chaplygin model in self-decelerating DGP gravity 
have been studied elsewhere (Bouhmadi-Lopez & Lazkoz l2007l l. 
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Fig. 1. The closed contour is the confidence region in the (Qj],, Q,.)- 
plane from a fit to SNela data. The point of best fit is marked by a 
cross. The curves crossing from left to right correspond to the flat- 
space condition Eq. (Hi for the upper lcr value (top), the central val- 
ue (middle) and the lower lcr value (bottom) of the parameter Q^. 
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